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Abstract 

We analyze the B — ► (J/tp, XcO, Xci, n c )K^ decays, which have small or even vanishing 
branching ratios in naive factorization assumption (FA) . We calculate nonfactorizable cor- 
rections to FA in the perturbative QCD approach based on Ut factorization theorem. The 
charmonium distribution amplitudes are inferred from the non-relativistic heavy quarko- 
nium wave functions. It is found that the nonfactorizable contributions enhance the 
branching ratios and generate the relative phases among helicity amplitudes of the above 
modes. Most of the observed branching ratios, polarization fractions, and relative phases, 
except those of B — > r/ c K, are explained. Our predictions for the B — > (xco, Xci, Vc)K* 
decays can be compared with future data. 



* chchen@mail . ncku . edu . tw 
thnli@phys. sinica. edu.tw 



1 INTRODUCTION 



It has been known that the naive factorization assumption (FA) [1, 2] does not apply to exclusive 
B meson decays into charmonia, such as B — > J/ipK [3]. These modes belong to the color- 
suppressed category [4], for which predictions from FA are always small due to the vanishing 
Wilson coefficient a 2 ~ 0. However, the branching ratios measured by Babar recently [5] 

B(B + ^J/^jK + ) = (10.61 ± 0.15 ± 0.48) x KT 4 , 
B(B° -> J/i/jK ) = (8.69 ± 0.22 ± 0.30) x 10" 4 , (1) 

imply a larger parameter a,2(J/ipK) ~ 0.20 — 0.30 [3]. The same difficulty has appeared in other 
similar decays B — > (xco? Xci, Vc)K: the observed branching ratios are usually many times larger 
than the expectations from FA. For example, the B — > XcoK decays do not receive factorizable 
contributions, but the data of their branching ratios are comparable to those of B — > J/ipK in 
Eq. (1). 

Many attempts to resolve this puzzle have been made in more sophisticated approaches 
(for a review, see [6]). The large a,2(J/ipK) leads to a natural conjecture that nonfactoriz- 
able contributions, such as the vertex and spectator corrections from Figs. l(a)-l(d) and from 
Figs. 1(e) and 1(f), respectively, may play an important role. It has been found [7] in the QCD- 
improved factorization (QCDF) [8] that the above nonfactorizable contributions, resulting in 
the branching ratio B(B° — ► J/ipK°) « 1 x 10~ 4 , are too small to explain the data. Since only 
the leading-twist (twist-2) kaon distribution amplitude was included in [7], the authors of [9] 
added the twist-3 contribution, which is chirally enhanced, though being power-suppressed in 
the heavy quark limit. Unfortunately, Figs. 1(e) and 1(f) generate logarithmical divergences 
from the end-point region, where the parton momentum fraction of the kaon is small. To make 
an estimation, arbitrary cutoffs for parameterizing the divergences, i.e., large theoretical uncer- 
tainties, have been introduced. The end-point singularities become more serious in the QCDF 
analysis of the B — > (xco, Xci)K decays [10]. It has been claimed that the data of those modes 
involving J/ip, Xci, and i] c are not understandable in QCDF, and only those involving Xco are 
[7,10,11]. 

The B — > J/ipK decays have been also studied in light-cone sum rules (LCSR) [12, 13], and 
the small branching ratio from the factorizable contribution [12], 

B(B -> J/^K) iact « 3 x 10^ 4 , (2) 

was confirmed. The nonfactorizable contributions considered in [12, 13] refer, however, to those 
arising from the three-parton kaon distribution amplitudes. In the corresponding diagrams an 
additional valence gluon from the kaon attaches one of the charm quarks in the J ftp meson. 
Adding these pieces, one arrived at the parameter a2{J/ipK) ~ 0.14 — 0.17, which is still 
insufficient to account for the data. The calculation of Fig. 1 in LCSR, involving two-loop 
integrals, has not yet been performed. Hence, the conclusion in [12, 13] simply indicates that 
the nonfactorizable contributions to the B — > J/ipK decays from higher Fock states of the 
involved mesons are negligible. 

We doubt that the conclusion from QCDF is solid for two reasons. First, the end-point sin- 
gularities render the estimation of the nonfactorizable contributions out of control. Second, the 
simple asymptotic models were adopted for the charmonium distribution amplitudes, which, 
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Figure 1: Feynman diagrams for nonfactorizable corrections to B — > J/ipK. 

without any theoretical and experimental bases, may underestimate the nonfactorizable contri- 
butions. Motivated by the above observations, we shall investigate the nonfactorizable correc- 
tions to FA from Fig. 1 using the following method. We take into account the vertex corrections 
through the variation of the renormalization scale for the factorizable contributions, because 
their evaluation requires more theoretical inputs, which have not yet been available. The spec- 
tator amplitudes will be computed in the perturbative QCD (PQCD) approach based on kx 
factorization theorem [14, 15, 16, 17]. The formalism, similar to that of the color-suppressed 
mode B° — > D 7t° [18], is free from the end-point singularities. We infer the charmonium dis- 
tribution amplitudes from the non-relativistic heavy quarkonium bound-state wave functions, 
which have been shown to well explain the cross section of charmonium production in e + e~ 
collisions [19]. Based on the universality of hadron distribution amplitudes, these more sophis- 
ticated models should be able to account for the exclusive B meson decays into charmonia, if 
QCD factorization theorem works. 

The B — > J/ipK^ decays have been studied in PQCD [20, 21, 22], compared to which 
the new ingredients of this work are: 1) we clarify several controversial statements on the 
end-point singularities in the B — > (J/ip,x c i)K analysis in the literature; 2) the meson 
distribution amplitudes derived from QCD sum rules [23, 24] are included up to two-parton and 
twist-3 level; 3) the models for the charmonium distribution amplitudes have both theoretical 
and experimental bases; 4) not only the B — > J/ipK^*' decays, but B — > (xco,Xci,Vc)K^ 
are investigated. It will be shown that the obtained branching ratios, polarization fractions, 
and relative phases among various helicity amplitudes, except those of B — > rj c K, are all in 
consistency with the existing data. The observed B — > r] c K branching ratios, significantly 
larger than our results, are the only puzzle. Therefore, our conclusion differs from that drawn 
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in QCDF [7, 10, 11]. The predictions for the B — > (xdo, Xci, Vc)K* modes can be compared with 
future measurements. 

We present our formalism for the B — > J/ipK^ decays as an example in Sec. II. The 
factorizable contributions are treated in FA, since the B — > transition form factors in 
the present case are characterized by a low scale, and may not be calculable. The vertex and 
spectator corrections are handled in the way stated above. The same formalism is then applied 
to the other exclusive B meson decays into charmonia in Sec. III. Section IV is the conclusion. 
The Appendix collects the expressions of the involved meson distribution amplitudes and of 
the factorization formulas for all the spectator amplitudes. For a review of the studies of 
semi-inclusive B meson decays into charmonia, refer to [25]. 



2 THE B J/tpK^ DECAYS 

We write the B (J/ip) meson momentum Pi (P 2 ) in the light-cone coordinates as 



P 1 = ^|(1,1,0 T ), P 2 = ^|(l,r 2 2 ,0 T ) , (3) 

with the mass ratio r 2 = mj^/m^. The meson momentum is then given by P 3 = Pi — P 2 . 
The polarization vectors of the J/ip meson are parameterized as 

e 2L = -jL- (l, -r 2 2 , T ) , e 2T = (0, 0, 1 T ) . (4) 
The relevant effective Hamiltonian for the B — > J/ipK^ decays is 

HcS = 7§ y* v cs\ c MOi + cmo 2 ] - v tb v* c k bi)o k j , (5) 

with the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V and the four-fermion opera- 
tors, 



Oi = 


(SiCj 


)v-A(Cjbi) V -A , 2 = (SiCi) V -A(Cjbj)v-A , 


o 3 = 




)v-A J2{qjQj)v-A , 4 = {sibj) V - A Y(qjqi)v-A , 
q q 


o 5 = 




)v-A J2(QjQj)v+A , Oq = (Sib^v-A J2(QjQi)v+A , 

q q 


o 7 = 




3 

h)v-A e ?fe<?j)v+A , 8 = -(sibj) V -A e q(9jQi)v+A ■ 
q A q 


o 9 = 




3 

b^v-A^eqfaq^v-A , O w = -(sib^v-A^egfaq^v-A 

q 1 q 



(6) 

i, j being the color indices. Below we shall neglect the tiny product V u bV* s , and assume 
VcbV* s = —V t bV t * s for convenience. 
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2.1 Factorization Formulas 



The B — > J/ipK decay rate has the expression, 

1 



32tt 



^|K b | 2 |y cs | 2 m 3 B (l-r 2 2 )Vl- (7) 



Computing the factorizable contribution to the amplitude A in PQCD, we found that its 
characteristic scale is around 1 GeV. With such a small hard scale, the perturbation theory 
may not be reliable. Therefore, we do not attempt to calculate the factorizable contribution, 
but parameterize it in FA. The amplitude A is then written as 

A = ^/j/^K^+M^, (8) 

with the J ftp meson decay constant fj/^ and the spectator amplitude MS J I^ K>I . The form 
factor Fi(q 2 ) is defined via the matrix element, 



{K{Pz)\h»s\ B ( P i)) = F ^ 2 ) 



m B ~ m K 
IT 



(Pi + " B q2 K % 



,2 



q = P\ — P3 being the momentum transfer, and itlk the kaon mass. The effective Wilson 
coefficient a e ff sums the contributions from both the tree and penguin operators in Eq. (5): 

a eS (n) = a 2 (fj) + a 3 (fx) + a 5 (n) , (10) 

with 

C 2 



a 2 = Ci + 



C4 . 3 / Ciq 



„ C4 3 /„ 



and the charm quark charge e c = 2/3. 

After parameterizing the factorizable contribution in FA, the soft dynamics is absorbed into 
the form factor F\. The remaining piece, i.e., the Wilson coefficient, is dominated by hard 
dynamics of 0(m b ), m b being the b quark mass. We shall set the renormalization scale to 
ji — m b 4.4 GeV, which corresponds to a e s(m b ) « 0.1. The spectator amplitude M-VM is 



characterized by a scale of 0(y Am b ) with A being a hadronic scale [18, 26, 27]: the hard gluons 
in Figs. 1(e) and 1(f), carrying the difference between the momenta of the soft spectator in the 

B meson and of the energetic spectator in the kaon, are off-shell by 0(\jKm b ). Strictly speak- 
ing, once a B meson transition form factor is treated as a soft object, the vertex corrections in 
Figs. 1(a)- 1(d) should be added, which become of the same order as the spectator amplitude. 
This is the counting rules of QCDF. The PQCD formalism for the vertex corrections requires 
the inclusion of the transverse momenta kr of the charm quarks, because of the end-point 
singularities in some modes. The kr- dependent charmonium wave functions are then the nec- 
essary inputs, which, however, have been determined neither theoretically nor experimentally. 
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Figure 2: Dependence of a c g on the renormalization scale //. The unit for the horizontal axis 
is GeV. The solid, dashed and dash-dotted lines stand for a eff without the vertex corrections, 
the real part of a e g- with the vertex corrections, and the imaginary part of a e ^ with the vertex 
corrections, respectively. 



Therefore, the vertex corrections appear as a theoretical uncertainty eventually, which can be 
covered by a variation of the scale \x between 0.5mfc and 1.5m.& as shown below. 

In terms of the notation in [28] , we decompose the nonlocal matrix elements associated with 
longitudinally and transversely polarized J/ip mesons into 

(JMP,e L )|c(z) iC (0),|0> = f 1 dxe" Pz \m J/ ^ L ] l ^ L (x) + y L P^x) 



(J/^(P,e T )|c(z) jC (0) 2 |0) 



V2Nc 
1 



£dxe ixP - z {mj /f yTh* V (x) + \fa W 



(12) 
(13) 



\f2~W c jo 

respectively, which define the twist-2 distribution amplitudes \I/ L and \I/ T , and the twist-3 
distribution amplitudes \&* and with the c quark carrying the fractional momentum xP. 
It is confusing that both the structures fa and fa f* were regarded as being leading-twist in 
the analysis of the B — ► J/ipK decays [7, 9]. In fact, it is fa P that is leading- twist, which 
contributes only to the B — > J/ ipK* decays. 

The vertex corrections to the B — ► J/ipK decays, denoted as fj in QCDF, have been 
calculated in the NDR scheme [7, 9]. Their effect can be combined into the Wilson coefficients 
associated with the factorizable contributions: 



a 2 



a 3 



a 5 



C + — + 
1 N r 4tt N, 



a s C F 

L-2 



,-, , a s C F 

( ' n ■ — + Tk W c Ca 
■Or 



9± 



_18 + I21n^ + /j 
-18 + 121n^ + // 



^5 + ^4 



6 - 12 In 



fi 



with the function. 

Si 



2^/2N~ c 



fj/i> 



dx 3 ty L {x 2 ) 



3(1 - 2x 2 
l—x 2 



In 22 — 37n + 3 ln(l — r\) + 



2r 2 2 (l-x 2 ) 
1 — r\x 2 



(14) 



(15) 
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where those terms proportional to r\ have been neglected. We have also neglected the contri- 
butions from the electroweak penguin operators for simplicity. As shown in Fig. 2, including 
the vertex corrections leads to a larger |a e ff(mb)| ~ 0.12, which reproduces Eq. (2) roughly. We 
have confirmed that the variation of \i between 0.5rrib and 1.5mb is wide enough for taking into 
account the effect of the vertex corrections. 

The definitions of the B meson wave function <3> B and of the kaon distribution amplitudes 
§ K are referred to [27, 28]. It is then straightforward to derive the spectator amplitude, 



where the amplitudes M.^^^ and M.^^^ 



(J/^K) 



+ M 



(J/ijjK) 



(16) 



result from the (V-A)(V-A) and (V-A)(V+A) 
operators in Eq. (5), respectively. Their factorization formulas are given by 



M VI*K) 



16nC F ^2N c j\dx] J™ b l db l <^ B (x l , h) 

x{ [(1 - 2r 2 2 )(l - x 2 )^ K (x 3 )^ L (x 2 ) + Ki$ K (x 3 )¥(x 2 ) 

-r K (l - r 2 2 )x 3 $ p K (x 3 )^ L (x 2 ) + r K (2rf (1 - x 2 ) + (1 - r 2 2 )x 3 ) ^ K (x 3 )^ L (x 2 ) 

xE A {tf)hf{x 1 ,x 2 ,x 3 ,b 1 ) 

- \x 2 + (1 - 2r 2 2 )x 3 )$ K (x 3 )V L (x 2 ) + r 2 2 {2r K ^ K {x 3 ) - ^k(x 3 ))¥(x 2 ) 

-r K (l - r 2 2 )x 3 <S> p K (x 3 )ty L (x 2 ) - r K {2r\x 2 + (1 - r 2 2 )x 3 ) ^ K (x 3 )^ L (x 2 ) 

xE A {tf)hf (x^x 2 ,x 3 M) , 
16nC F ^2N c f\dx] b 1 db 1 ^ B {x 1 , h) 



(17) 



x 



1 - x 2 + (1 - 2rl)x 3 )<b K (x 3 )m L (x 2 ) + r 2 2 (2r K ^ K (x 3 ) - -<f> K (x 3 ))*\x 2 ) 



{ 

-r K (l - r 2 2 )x 3 <5> p K (x 3 )V L (x 2 ) - r K (2r 2 2 (l - x 2 ) + (1 - r 2 2 )x 3 ) <f>° K (x 3 )V L (x 2 ) 
xE 6 (t { J ) )h^(x 1 ,x 2 ,x 3 ,b 1 ) 

- (1 - 2r 2 2 )x 2 <5> K {x 3 )q L {x 2 ) + K- 2 2 <5> K (x 3 )m\x 2 ) 

-r K (l - r 2 2 )x 3 & p K (x 3 )m L (x 2 )+r K [2r\x 2 + (1 - r\)x 3 ) §° K {x 3 )^ L \x 2 ) 
xE & {tf)hf{ Xl) x 2) x 3 M)} , (18) 

with the color factor C F = 4/3, the number of colors N c = 3, the symbol [dx] = dx\dx 2 dx 3 and 
the mass ratio r K = uIq /ttib, being the chiral scale associated with the kaon. 
The evolution factors are written as 



Ei(t) = a s (t)a[(t)S(t)\ 

63=61 ) 



(19) 



with the Wilson coefficients, 



C2 2. 



C 4 + -e c C 10 



a R = 



N r 



C 6 + -e c Ct 



(20) 



The Sudakov exponent is given by 



S(t) = S B (t)+S K (t) 
S B (t) = exp 



s(xiP 1 + ,6i) - \ [ — 7(a s (^)) 

3 Jl/bi pL 



s K (t) 



exp 



-s(x 3 P 3 -, b 3 ) - s((l - x 3 )P 3 ~, b 3 ) - 2 f ^Mti) 



(21) 
(22) 



with the quark anomalous dimension 7 = —a s /ir. Note that the coefficient 5/3 of the quark 
anomalous dimension in Eq. (21) differs from 2 in Eq. (22). The reason is that the rescaled 
heavy-quark field adopted in the definition of the heavy-meson wave function has a self-energy 
correction different from that of the full heavy-quark field [29]. The variables b\ and 63, con- 
jugate to the parton transverse momenta hyp and k 3 T, represent the transverse extents of the 
B and K mesons, respectively. The explicit expression of the exponent s can be found in 
[30, 31, 32]. The above Sudakov exponentials decrease fast in the large b region, such that the 
B — > J/ r ipK spectator amplitudes remain sufficiently perturbative in the end-point region of 
the momentum fractions. 

The hard functions j = 1 and 2, are 



h 



K (^~D~jm B bi) - K (y/Dm B bi) 
D-Dj\ fH^ l^D~\m B b 1 ) - KoiVDmsh) 

with the arguments, 

D = xix 3 (l -rl), 

Di = (1 - x 2 )xirl + (x! + x 2 - l)x 3 (l - rl) + -- 

D 2 = xix 2 rl + {xi - £2)2:3(1 - r 2 ) + Q - xfj rl . 

The hard scales t are chosen as 

= max(\^Dm B , ^\Dj\m B , l/&i) . 



for Dj > 
for Dj < 



x 2 f 



2 1 



(23) 

(24) 
(25) 

(26) 
(27) 



Without the transverse momenta kr, the terms containing 2r^r|$^ : 1 ]/* in Eqs. (17) and (18) 
are logarithmically divergent due to the end-point singularities from x 3 — > [9] . Since the end- 
point singularities appear only at the power of r K r\, instead of at leading power, the spectator 
contributions are expected to be evaluated more reliably than the factorizable contributions in 
PQCD. 

For the B — > J/ipK* modes, the parametrization of the kinematic variables in Eqs. (3) 
and (4) still apply. We further define the polarization vectors e 3 for the K* meson, which is 
orthogonal to the K* meson momentum P 3 . In this case the K* meson mass should be kept in 
the intermediate stage of the calculation, and then dropped at the end according to the power 
counting rules in [27]. The decay rate is given by 



32tt 



\ 2 ml 



;) 3 £-4 (CT) l4 (CT) . 



(28) 
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The amplitudes for different final helicity states are expressed as 



A (a) = -\a e Mfj /lP F(m* )+M 



(J/^K*) 



e 2T ' e 3T 



2 "+M 



■a/3-yp* * 
1 e e 2a e 3(3 



r 2 a e fi(^)fj/^A l (mj / ^ 



N 



m. 



r 2 a e s(n)fj/ 1 pV(m J/iJ ) + M> T 



with the combination of the form factors, 

l + r K 



(29) 



(30) 



2r K * " _iV " Jlip/ 2r K *{l + r K *)~ 

and the mass ratio r K * = m K */m B . The first term in Eq. (29) corresponds to the configuration 
with both the vector mesons being longitudinally polarized, and the second (third) term to 
the two configurations with both the vector mesons being transversely polarized in the parallel 
(perpendicular) directions. 

The effective Wilson coefficient a e ff is the same as in Eq. (10). The B — > K* transition form 
factors Ai, A 2 , and V are defined via the matrix elements, 

2iV(q^ 



(/r(P s ,e 3 )|67 |, a|S(Pi)> 
(ir(P 3 ,e3)|^75s|50Pi)> 



; e Hv^zp^ic 

tub + rriK* 

2m K ,A ( y q 2 )^-q" + (m B + m^A^q 2 ) 



(31) 



-A 2 (q 2 ) 



£3 



niB + rriK* 



Pi + P3 ~ 



m. 



m 



Q 



The spectator amplitudes M.^n,t ^ from Figs. 1(e) and 1(f) are written as 

vl LA,NA,TA "r JVI L6,N6,T6 > 



JVI L,N,T 



(32) 



(33) 



where the explicit factorization formulas for the amplitudes M^l1^n1,ta an d for M^ll^m^e, are 
presented in the Appendix. 

The helicity amplitudes are then defined as, 

'aMfj/^Fim^+M^^ 



A L 

A, 



G 



-GV2 r^Mfj/^A^m^+M^ 



(J/tpK*) 



(34) 



A ± = -Gr 2 r K ^2[{v 2 ■ v z f - 1] [r 2 a cff ( / i)/ J/ ^(m 2 w ) + M^ K ' ] ] , 

for the longitudinal, parallel, and perpendicular polarizations, respectively, with the velocities 
v 2 = P 2 /mji^ and 1*3 = P 3 /mx*. The normalization factor G has been chosen such that the 
following relation is satisfied: 

\A L \ 2 + \A\\\ 2 + \Ax\ 2 = I . (35) 

In our convention the relative phase 5\\ (5±) between A\\ (A±) and A L takes the value of n in the 
heavy-quark limit, ie., as a cS is real and M^nt ^ vanish. We define the polarization fractions 
fh = \Al\ 2 , f\\ = \A\\\ 2 , and f± = \A_i_\ 2 . The B — > J/ipK* polarizations, significantly affected 
by the nonfactorizable contributions, were not discussed in [33]. 
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2.2 Numerical Analysis 

To perform the numerical analysis, we need the information of the involved meson distribution 
amplitudes. The B meson wave function <3> s and the K^> meson distribution amplitudes have 
been studied intensively, whose expressions are collected in the Appendix. The charmonium 
distribution amplitudes attract less attention in the literature. In [34] they were assumed to be 
identical to the corresponding light-meson distribution amplitudes: those of the J / if) (r] c ) meson 
are the same as of the p (tt) meson. Then the twist-3 charmonium distribution amplitudes do 
not vanish at the end points x — 0, 1, where the valence charm quarks become highly off- 
shell. Such a functional shape seems to contradict to the intuition. More realistic asymptotic 
models have been proposed recently in [19] , which were inferred from the non-relativistic heavy 
quarkonium bound-state wave functions with the same quantum numbers. We shall adopt 
these J /if) and r\ c meson distribution amplitudes, and derive the Xco and Xci meson distribution 
amplitudes following the similar procedure in the next section. 

The J /if) meson asymptotic distribution amplitudes are given by [19] 



¥(x) 

V v (x) 



* r ( x ) = 9.58-^=^(1 - x) 



10.94 
1.67 



2y/2N c 

fj/ ^ (l-2x) 2 



x(l — x) 



1 - 2.&r(l - x) 



0.7 



x(l — x) 



0.7 



1 + (2x - 1) 



1 - 2.8x(l -x 

x(l — x) 



l-2M(l-x) 



0.7 



(36) 



in which the twist-3 ones \E r *' y vanish at the end points due to the additional factor [x(l —x)] ' 7 . 
Compared to [20] , we have distinguished the distribution amplitudes associated with the longi- 
tudinally and transversely polarized J /if} mesons, which exhibit the different asymptotic behav- 
iors. We do not distinguish the decay constant fj/^ for the normalization of \I/ L and fj^ for the 
normalization of The distribution amplitudes associated with the structures I (the identity) 
and 7^75 (the pseudo- vector) diminish like 1 — 2m c /mj/^, m c being the c quark mass. The 
above models have been shown to yield the observed cross section of charmonium production 
in e + e~ collisions [19]. If factorization theorem works, the same distribution amplitudes, due 
to their universality, should be able to explain the exclusive B meson decays into charmonia. 





F(0) a b 




0.35 1.58 0.68 


F (q 2 ) 


0.35 0.71 0.04 


V(q 2 ) 


0.31 1.79 1.18 


A (q 2 ) 


0.31 1.68 1.08 


A^q 2 ) 


0.26 0.93 0.19 


A 2 (q 2 ) 


0.24 1.63 0.98 



Table 1: Parameters for the B — > transition form factors. 
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X(J PC ) 


m x [GeV][35] 


/x[MeV] 




2.980 


420 ± 50 [36] 


JM1-) 


3.097 


405 ± 14 [12] 


Xco(0 ++ ) 


3.415 


360 [37] 


Xcl(l ++ ) 


3.511 


335 [37] 



Table 2: Decay constants and masses of various charmonium states. 

For the B — > transition form factors, we employ the models derived from the light-front 
QCD [38], which have been parameterized as 

F(g2) = 1 - a(q 2 /m 2 B ) + b(q 2 /m 2 B ) 2 ' (37) 

with the constants F(0), a, and b being listed in Table 1. Using the inputs in Table 2 [39] for 
the various charmonium states, the CKM matrix elements V c b = 0.040 and V cs = 0.996, the 
quark mass m t = 174.3 GeV for the Wilson coefficients, the lifetimes t b o = 1.56 x 10~ 12 sec 
and t b ± = 1.67 x 10~ 12 sec, and the Fermi constant G F = 1.16639 x 10~ 5 GeV -2 , we derive 
the branching ratios, 



B{B + - 


-> j/i)K + ) -- 


- (9.20i 6 7 ;^) x 10- 4 


B(B° - 




(8.60^) x 10" 4 , 


B{B + - 


-> j/^k* + ) 


= (9-95^; 2 6 ) x 10 -4 


B(B° - 


-> J/^K*°) = 


= (9.30+ 4 ;^) x 10- 4 



(38) 

We have also derived the the polarization fractions and the relative phases among the helicity 
amplitudes for the B — > J/ipK* decays, 

A = o.73±8:Si, /n = o.i7^, /x = o.io±g:S , 

5|l=2.57l^ 2 , tfj. = 2.53±g:Sg - (39) 

Both the B + — > J/ipK* 0+ and B° — > J/if)K*° modes have the same polarization fractions, since 
they differ only in the lifetimes in our formalism. The theoretical uncertainties come from the 
variation of the renormalization scale \x. We have checked the sensitivity of our results to the 
choice of the b quark mass: m b = 4.8 GeV increases the branching ratio B(B° — > J/ipK°) only 
by 10%. 

The recent Babar measurement gave [5] 

B{B + -> J/ipK* + ) = (14.54 ± 0.47 ± 0.97) x 10~ 4 , 

B(B° -> J/ipK* ) = (13.09 ± 0.26 ± 0.77) x 10~ 4 (40) 

and the Belle measurement gave [40] 

f L = 0.585 ± 0.012 ± 0.009, /y = 0.233 ± 0.013 ± 0.008, f ± = 0.181 ± 0.012 ± 0.008 , 
6\\ = 2.888 ±0.090 ±0.008 , 5± = 2.903 ± 0.064 ± 0.010 . (41) 
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It is found that the consistency between the theoretical and experimental values for the B — > 
J/ipK^ branching ratios is reasonable. The predicted longitudinal polarization fraction for 
the B — > J/ipK* is a bit larger, and the predicted relative phases are a bit smaller than the 
data. However, we point out that the polarization fractions and the relative phases could be 
modified by adding the higher Gegenbauer terms to the J / ip meson distribution amplitudes, and 
that the consistency will be improved. This fine tuning will not be performed here. Roughly 
speaking, it is very promising to understand the B — > J/ipK^ data in the PQCD framework. 

Another remark is as follows. Most of the J/ip meson distribution amplitudes, except 
exhibit maxima at x = 1/2, such that the valence charm quark, carrying the invariant mass 
X let ~ m 2 , is almost on-shell. However, \1>* has two humps with a dip at x = 1/2. We argue 
that the models in Eq. (36) are the asymptotic ones. Including the higher Gegenbauer terms, 
such as those for a light vector meson, 



¥(x) = 10.94-^={(l-2:r) 2 + a^(2:r - l) 2 [5(2:r - l) 2 -3] 



2^/2N~ c 

+a J 2 ^[3 - 30(2x - l) 2 + 35(2x - l) 4 ]} 



/, x -10.70 

x(l — X) 



1 - 2.8x(l - x) 



(42) 



possess three humps with the major one located at x = 1/2. Adopting this parametrization, 
the Gegenbauer coefficients and aj^ should be determined from the data of charmonium 
production in e + e~ annihilation, before it is applied to the exclusive B meson decays into 
charmonia. 



3 B - (xco, Xci, Vc)KU DECAYS 



Before calculating the spectator contributions to other similar modes, we clarify a confusing 
statement in the literature. It has been claimed [10] that the end-point singularities from 
x 3 — > cancel between the spectator diagrams Figs. 1(e) and 1(f) in the QCDF formalism for 
the B — > J/tpK decays, but do not for the B — > x c \K decays. Our opinion is that the end-point 
singularities should appear in both modes at the power of r|, if no approximation is made. The 
singularities cancel in [9] because of the relation, 



m. 



2x\ , (43) 



where the decay constant fj,^ is associated with the normalization of the twist-3 distribution 
amplitude This relation arises from the exact on-shell condition of the valence charm 
quarks. The singularities were present in [10], since the authors, considering only the leading- 
twist Xci meson distribution amplitude, got no chance to apply Eq. (43). We agree on the 
comment made in [9]: it is not clear why the singularities cancel in [7] for the B — > J/ifjK 
decays, viewing that they treated the ratio on the left-hand side of Eq. (43) as a constant. The 
exact on-shell condition is indeed not necessary, because a deviation from mass shell, as long as 
being power-suppressed, is allowed. Furthermore, the momentum fraction x 2 , running between 
and 1, is hardly related to the left-hand side. Once Eq. (43) is not postulated, the end-point 
singularities exist in the terms proportional to r\ for the B — > J/ipK decays [20], and also for 
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the B — > Xc\K decays. Note that these singularities disappear in the PQCD approach based 
on kx factorization theorem. 

The PQCD analysis of the B — > XciK^ decays is similar to that of B — > J /ipK^*\ To cal- 
culate the nonfactorizable spectator amplitudes, we consider the Xci meson distribution ampli- 
tudes defined via the nonlocal matrix elements associated with the longitudinal and transverse 
polarizations, 

(Xci(P,CL)|c(«) iC (0),|0> = -^J^ dxe ixPz \m Xcl [ lb ^}ijXi( x ) + [75 h ? ]yxi(z)},(44) 

( Xcl (P,e T )|c(^c(0) z |0> = -J=£dxe^{m x J l5 MnxX{ x ) + [75 fr f]^(x)}(45) 

respectively, with the Xci meson mass m Xcl . The asymptotic models for the twist- 2 distribution 
amplitudes Xi( x ) an d Xi( x )i an d for the twist-3 distribution amplitudes Xi( x ) an d Xi ( x ) wn l 
be derived following the prescription in [19]: they can be extracted from the n — 2, I — 1 
Schrodinger states for a Coulomb potential, whose radial dependence is given by 

*Sch(r) oc rexp(-g B r) , (46) 

qs being the Bohr momentum. Fourier transformation of the above solution leads to 

P - la 2 

*»<*> K (FT# • (47) 



with k 2 = |k| 2 . Employing the substitution [19], 

k, - k, . /.•: --(.'• -r!^. ,,,7, "' • (18) 



_ x m jD 2 m 2 + fcy 



x (x = 1 — x) being the c (c) quark momentum fraction, we obtain the heavy quarkonium 
distribution amplitude, 

«(x) ~ / Jfa-W*. fa.) oc xx { gli=g(i±j I , (49 ) 
with f 2 = q 2 B /m 2 c . 

To impose a fractional power 1 — v 2 to the factor in the curved brackets [19], we neglect the 
v 2 term in the numerator. Otherwise, the numerator may become negative for 



1 V 1 ~ T+b 5 <2x<1 + \j 1 - ' (50) 

and the analyticity of the distribution amplitude will be lost. We then propose the Xci meson 
distribution amplitudes inferred from Eq. (49), 

XlW oc^ W { Tr ^g^_ }'", (51) 

with $ asy (x) being set to the asymptotic models of the corresponding twists for light vector 
mesons. It is now clear that after imposing the fractional power 1 — v 2 , Eq. (51) approaches 
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the light meson distribution amplitudes as v 2 — > 1, and the heavy quarkonium distribution 
amplitudes the same as Eq. (49) as v 2 — > [19]. Fixing the normalization, we derive 

f f \N A /-, M A 0.7 

JXcl 



X{ (x) = 27.46 

fxci 



15.17 



3.60 



fxci 



2^2N, 
1 - 2x) 2 



x(l — x) 



x(l -x)[l -4x(l -x)] 
[1 - 2.8x{l-x)] 2 



x(l -x)[l -4x(l -x)} 



0.7 



1 + (2x- if 



[1 -2.8x(l -x)} 2 j ' 

f x(l -x)[l -4x(l -x)] 
\ [1 - 2.8x{l - x)} 2 



0.7 



(52) 



where the same decay constant f Xcl has been assumed for the longitudinally and transversely 
polarized Xci meson. Note that the dip at x = 1/2 is a consequence of the P-wave Schrodinger 
wave functions. 

The B — > Xc\K^ decay amplitudes are written as 

A = aMf Xcl Fi(m 2 Xcl )+M^ K \ (53) 

,2 \ , kAXciK*) 



= -\aMf Xcl F{m 2 Xci )+Mf 



Z 2T ' t 3T 



■ afiip * * 
1 fc fc 2a fc 3/3 



m 



B 



r2f Xcl V(m x J+M$ 



(XciK*) 



(54) 



with the mass ratio r 2 



(55) 



m Xci/ m B> and the effective Wilson coefficient, 

a ef£ (n) = a 2 (n) + a 3 (/i) - a 5 (/i) . 

Similarly, we shall vary the renormalization scale \x in the range between 0.5m;, and 1.5m;,, 
which covers the effect of the vertex corrections. As in the case of B — > J/ipK, the infrared 
divergences in the vertex corrections cancel in the B — > x&K decays: the hard kernel of the 
vertex corrections is odd for the structure 75 jk under the exchange of x and 1 — x, while 
the corresponding distribution amplitude is even. The hard kernel vanishes for the structure 
75 / p. Express the spectator amplitudes M^ 01 ^ and M.^n% ^ as 

M (XdK) = M (XdK) + M (XdK) ? 

(56) 



M (XdK*) 
JVI L,N,T 



M (XdK*) , M (XdK*) 
- /w L4,Af4,T4 + JVI 1 



l L6,N6,T6 ' 



where the factorization formulas of the amplitudes M^ clK \ M^ 1 ^, M^na,t4 an d M^m,T6 



are shown in the Appendix. Their expressions are similar to those of B —. 
terms flipping signs. 

Using the inputs listed in Table 2 [39], we get 

B(B+ -> XdK + ) = (3.15±^D x 10- 4 , 
B(B° -> XciK ) = (2.9432) x 10" 4 , 
B(B+ -> XdK* 0+ ) = (2.99^) x lO" 3 , 
B(B° -> Xcl lT ) = (2.79l8;| 6 7 ) * 10~ 3 , 

/l = 0.381°:™ /I, = o.o7iHl, /± = o.55i8:8g, 



J/i/jK^ with some 



2.6H8:43, *± 



2.58 



+0.07 
-0.03 



(57) 
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The errors of the above predictions arise from the variation of the renormalization scale \i for 
the factorizable contribution. For the B — > x&K* branching ratios, the imaginary part of the 
spectator amplitude Ai^f clK dominates, such that the influence of the // dependence becomes 
mild. This is also the reason we obtain a large f±. The B — > XaK branching ratios in Eq. (57) 
are consistent with the recent Babar measurement [5], 

B(B + ^XdK + ) = (5.79 ± 0.26 ± 0.65) x 10" 4 , 
B(B°^XciK°) = (4.53 ± 0.41 ± 0.51) x 10~ 4 . (58) 

Note that these branching ratios, with a parametrization for the logarithmical end-point singu- 
larities, were estimated to be only about 10 -4 in QCDF [10]. Our predictions for the B — * Xc\K* 
decays, including the branching ratios, the polarization fractions (especially the large f±), and 
the relative phases, can be compared with data in the future. 

We then discuss the B — > XcoK decays, for which the nonfactorizable contributions dominate 
due to the absence of the factorizable ones. The nonlocal matrix element associated with the 
Xco meson is decomposed into, 

( Xc o(P)|c(z) jC (0);|0> = ^^^^{[fl^^+m^J/]^^)}, (59) 

which defines the twist-2 and twist-3 distribution amplitudes, Xo( x ) an d Xo( x ), respectively, 
m Xc0 being the Xco meson mass. To satisfy the identity (x c o(-P)|c7 /i c|0) = 0, Xo( x ) mus t be 
anti-symmetric under the exchange of x and 1 — x. Following the similar ansatz of constructing 
the Xci meson distribution amplitudes, we propose the following asymptotic models, 

f Xc0 . f x(l - x)[l - 4s(l - x)]\ - 7 



Xo(x) = 27.46 J ** (l-2x) 



I2s/W c y '{ [1 - 2.8x(l - x)\ 



2 



Xo[X) { 2V2N C \ [1 - 2.8x(l - x)f J ' [W) 

where the decay constant f Xc0 for the normalization of Xo an d Xo nas been assumed to be equal. 

In the QCDF analysis of the B — > (J/iJ;,Xci)K^ decays, the vertex corrections from 
Figs. l(a)-l(d) are infrared safe. However, for the B — > XcqK decays, the hard kernels of 
the vertex corrections are even for the structure /, and odd for f under the exchange of x and 
1 — x. The corresponding Xco meson distribution amplitudes are also even for the structure 
/, and odd for p. Therefore, the infrared divergences do not cancel [10]. To regularize these 
divergences, a binding energy 2m c — m Xc0 and a gluon mass m g have been introduced in [41] 
and [42], respectively. The vertex corrections to the B — > XcoK decays can also be handled in 
the PQCD approach. For the reasons stated in the Introduction, we do not attempt such a 
calculation here. Instead, we shall demonstrate that the spectator contributions are sufficient 
to account for the observed B — > XcoK branching ratios. The decay amplitudes are written as 

A = M(*cO*(*>) ^ M (XcoK(*)) = M (XcoK(*)) + M (X C0 KM) ^ (61) 

with the explicit factorization formulas of the amplitudes M^° K( ^ being referred to the Ap- 
pendix. There are also logarithmic end-point singularities in these amplitudes, if using QCDF, 
indicating that the PQCD approach based on kr factorization theorem is more appropriate. 
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Adopting the inputs in Table 2, we obtain the branching ratios, 



B{B+ 
B(B° 
B(B + 
B(B° 



■ XcoK + ) 
XcoK ) 



5.61 x 1(T 4 , 



-4 



5.24 x 10 
= 8.69 x 1(T 4 
= 8.12 x 10~ 4 . 



(62) 



There is no theoretical uncertainty from the variation of /i, since the factorizable contributions 
vanish in this case. Babar gave the upper bounds (the values) [43], 



and Belle gave [44] 



B(B + -> XcoK + ) < 8.9 (= 4.4 ± 3.3 ± 0.7) x 10~ 4 
B(B° -> XcoK ) < 12.4 (= 5.3 ± 5.0 ± 0.8) x 10~ 4 



B(B + -> XcoK + ) = (1.96 ± 0.35 ± 0.33^) x 10" 



(63) 



(64) 



where the third error comes from a model uncertainty. Hence, we conclude that the above data 
are understandable. The B + — > XcoK + branching ratio was found to be as small as 4.2 x 10" 5 
in [41] (with a large theoretical uncertainty from parameterizing the end-point singularities in 
QCDF) and (2 — 3) x 10~ 4 for the gluon mass m g = 0.5 — 0.2 GeV in [42]. Despite of several 
discrepancies in their factorization formulas for the vertex corrections and for the spectator 
amplitudes, the difference in their numerical results is not really essential due to the tunable 
parameters. 

The analysis of the B — * r) c K^ decays is similar to that of B — > XcqK^*\ The nonlocal 
matrix element associated with the Xco meson is decomposed into 



(ifc(P)|<fc);c(0),|0) 



dxe 



ixP-z 



[75 P]ijV v (x) +m Vc [-f 5 } lj r ] s (x) 



(65) 



which defines the twist- 2 and twist-3 r\ c meson distribution amplitudes, rf{x) and i] s (x), re- 
spectively, m V( , being the rj c meson mass. The asymptotic models for the r\ c meson distribution 
amplitudes have been determined in [19]: 



rf{x) = 9.58 
r] s (x) = 1.97 



fr)c 



2y/W c 



:X(1 — X) 



x(l — x) 



x(l — x) 



1 - 2.8x(l - x) 

0.7 



0.7 



1 - 2.8x(l - x) 



(66) 



where we have assumed the same the decay constant f Vc for the normalization of if and rf . 
The decay amplitudes are written as 



A = aMUF (rn 2 v ) + M^ K) , A = a^W^M^ + , 
for which the factorization formulas of the spectator amplitudes, 



(67) 



(68) 
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can be found in the Appendix. The B — > K transition form factor F and the B — > K* 
transition form factor A have been defined in Eqs. (9) and (32), respectively. The effective 
Wilson coefficient a e s(fi) for the B — > XcqK^ decays is also given by Eq. (55). Similarly, we 
shall vary the renormalization scale \x in the range between 0.5m;, and l.brrib, which covers the 
effect of the vertex corrections. The hard kernels of the vertex corrections are even for the 
structure 75, and odd for 75 f under the exchange of x and 1 — x. The corresponding r\ c meson 
distribution amplitudes are both even for 75 and 75 f. Therefore, the infrared divergences in 
the vertex corrections to the B — > r] c K decay, cancelling only for 75 f, still exist. 

Using the inputs in Table 2 [39] and the form factor parametrization in Eq. (37), we obtain 
the branching ratios, 



B{B + - 




) = (2.341^; 


) x 10" 4 


B(B° - 


- VcoK ) 


= (2.19l!;E) 


x 10~ 4 , 


B(B + - 


- VcK* + 


) = (2.8211:%; 


) x 10~ 4 


B(B° - 


- VcK*°) 


= (2-64±l:S) 


x 10- 4 . 



(69) 

Since the r\ c meson distribution amplitudes have no two humps, the B — > r\ c K branching ratios 
are expected to be smaller than the B — > J/ipK ones. However, the Belle measurement [45], 

B(B + ^r] c K + ) = (1.25 ± O.Ut° i° 2 ± 0.38) x 1(T 3 , 
B(B°^r] c K°) = (1.23 ± 0.23lg;^ ± 0.38) x 10" 3 , (70) 

and the Babar measurement [46], 

B(B + ^r] c K + ) = (1.34 ± 0.09 ± 0.13 ± 0.41) x 10" 3 , 
B(B°^r] c K°) = (1.18 ± 0.16 ± 0.13 ± 0.37) x 10" 3 , (71) 

are significantly larger than our predictions. The QCDF approach gave the B — * i] c K branching 
ratios about (1.4 — 1.9) x 10~ 4 for a rough estimate, which are also too small [11]. 

Because the r\ c meson distribution amplitudes in Eq. (66) have been shown to produce the 
observed cross section for charmonium production in e + e~ annihilation [19], they are supposed 
to explain the branching ratios of the exclusive B meson decays into charmonia according 
to the universality, even though these two processes involve dramatically different dynamics. 
Therefore, we conclude that the B — > rj c K decays are the only puzzle, and demand more studies. 



4 CONCLUSION 

The exclusive B meson decays into charmonia, B — > XK^*\ for X — J/ip, Xco, Xci, and r/ c have 
been analyzed in the QCDF approach in the literature. Among these modes, those with X = Xco 
were claimed to be explainable, and those with X = J/ifj,Xci, and i] c were not. However, 
this conclusion is not solid due to the serious end-point singularities in many QCDF decay 
amplitudes for the vertex corrections and for the spectator contributions, and due to the ad- 
hoc models of the charmonium distribution amplitudes. In this paper we have investigated these 
decays in an improved framework: the factorizable contributions are still treated in FA, since 
the involved B — > transition form factors, evaluated at the charmonium mass, may not be 
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calculated perturbatively; the effect of the vertex corrections was taken into account by varying 
the renormalization scale /i, because their estimation in the PQCD approach needs additional 
nonperturbative information, i.e., the distribution of the charm quark in its transverse degrees 
of freedom inside a charmonium; the spectator contributions were computed in PQCD, in 
which the end-point singularities are smeared by the Sudakov factor associated with the 
meson. The J/ip and r\ c meson distribution amplitudes have been inferred from the n — 1, 
/ = Schrodinger state for a Coulomb potential [19] , which produce the measured cross section 
of charmonium production in e + e~ collisions. To analyze the B — > (xco, Xci)K^ decays, we 
have obtained the Xco, Xci meson distribution amplitudes from the n — 2, I — 1 Schrodinger 
states following the similar procedure. That is, our models for the charmonium distribution 
amplitudes have been constrained theoretically and experimentally. 

Our investigation has indicated that only the B — > i] c K decays exhibit a puzzle, whose 
branching ratios are significantly smaller than the observed values. The data of the B — > 
(J/ip, Xco, Xa)K, J/ipK* decays, including the branching ratios, the polarization fractions, and 
the relative phases among the various helicity amplitudes, are all understandable within theo- 
retical uncertainty. If QCD factorization theorem works, the r\ c meson distribution amplitudes, 
explaining the charmonium production data, are supposed to explain the exclusive B meson 
decays into charmonia. Hence, the B — > i] c K modes require a more thorough study. We have 
predicted the branching ratios, the polarization fractions, and the relative phases associated 
with the B — > (Xci, Xco, Vc)K* decays, which can be compared with future measurements. 
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the National Science Council of R.O.C. under Grant No. NSC-93-2112-M-001-014. 



A FACTORIZATION FORMULAS 

For the B meson wave function, we adopt the model [16], 
$b(s, b) = N B x 2 {l - x) 2 exp 



4 


( xrriB s 


) 2 ^& 2 1 




V UJb y 





(72) 



with the shape parameter uj b = 0.4 GeV. The normalization constant N B = 91.784 GeV is 
related to the decay constant fs = 190 MeV (in the convention f w = 130 MeV). The K and 
K* meson distribution amplitudes have been derived from QCD sum rules [23, 24], 



& K (x 

$ K *(x 

& K *{x 



3/ 



K 



:X(1 



- x) {l + 0.51(1 - 2x) + 0.3[5(1 - 2xf - 1]} , 



= 1 J^= [l + 0.24C 2 1/2 (1 - 2x) - 0.11C 4 1/2 (1 - 2x) 
= (1 - 2s) [l + 0.35(10x 2 - 10s + 1) 

2y2N c l 



fl 



=x(l - s) [l + 0.57(1 - 2s) + 0.07C 2 3/2 (1 - 2s) 



(73) 
(74) 
(75) 
(76) 



— ^—(o.3(l - 2s) [3(1 - 2s) 2 + 10(1 - 2s) - ll + 1.68C 4 1/2 (1 - 2s) 
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0>^(s) 



-0.12s(l - s) + 0.36[1 - 6s - 2 ln(l - s)] J , 

Mi 



x(l -x)[l + 0.6(1 - 2s) + 0.04C 2 3/2 (1 - 2s) 

2V 2N c {l[ 1 + (! - 2a; ) 2 + O- 4 ^ 1 - 2a; ) 3 ] + 0.4C 2 1/2 (1 - 2s) 
+0.88C 4 1/2 (1 - 2s) + 0.48[2s + ln(l - x)} J , 
^={3(1 - 2s) [l + 0.19(1 - 2s) + 0.81(10s 2 - 10s + 1) 
-1.14s(l -s) +0.48[l-6s-21n(l-s)]} , 



+0.06(1 - 2s) 2 [5(1 - 2s) 2 - 3] + 0.36 [1 - 2(1 - 2s)(l + ln(l - s))] J , (77) 
3(1 - 2s) [l + 0.2(1 - 2s) + 0.6(10s 2 - 10s + 1)" 

(78) 
(79) 

(80) 
(81) 
(82) 



with the Gegenbauer polynomials 

Cl /2 (0 = ^e-l), Q 1/2 (O = ^(35£ 4 -30£ 2 + 3), C 2 3/2 (0 = ^(5£ 2 -l) 



and the decay constants fx = 160 MeV, fx* = 200 MeV, and f^, = 160 MeV. The coefficients 
of the Gegenbauer polynomials correspond to the masses ttik = 0.49 GeV and =1.7 GeV. 
We adopt the K* meson mass m K * = 0.89 GeV. 

A.l B^(J/*,Xci)K^ 

We present the spectator amplitudes for the B — > (J/ijj,Xci)K^ decays below, where the 
symbol X represents the charmonium J/ip or Xci- 



M X A K 



lQnC F \FlN c [\dx] H bxdbx^BixiM) 
v Jo Jo 

x{ [(1 - 2r 2 )(l - x 2 )<$> K {x 3 )X L {x 2 ) + C Y Kl^ K {x^)X t {x 2 ) 
-r K (l - r 2 2 )x 3 <5> K (x 3 )X L (x 2 ) + r K (2r 2 (l - s 2 ) + (1 - r 2 )s 3 ) <^(s 3 )X L (s 2 ) 
xE i {t { p)h { p{x l) x 2) x^b l ) 

- \x 2 + (1 - 2r 2 )s 3 )$ x (s 3 )X L (s 2 ) + r 2 2 (2r K $ K (x 3 ) - ^ K (x 3 ))X t (x 2 ) 
-r K (l - r 2 )s 3 $^(s 3 )X L (s 2 ) - r K (2r 2 s 2 + (1 - r 2 )s 3 ) <S> K (x 3 )X L (x 2 ) 
xE 4 (t^)hf\x 1 ,x 2 ,x 3 ,b 1 )} , (83) 



Mq K 



CxltoC F yf2N c £[dx] hdb&Bfa, 61) 
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x{ [(1 - x 2 + (1 - 2r 2 2 )x 3 )$ K (x 3 )X L (x 2 ) + ( x r 2 2 (2r K ^ K (x 3 ) - U> K (x 3 ))X\x 2 ) 
-r K (l - r 2 2 )x 3 $ p K (x 3 )X L (x 2 ) - r K (2r 2 2 (l - x 2 ) + (1 - r 2 2 )x 3 ) ^ K (x 3 )X L (x 2 ) 
*Es(t { p)h { p(xi ) x 2) x 3 M) 

- "(1 - 2r 2 2 )x 2 <b K {x 3 )X L (x 2 ) + K 2 2 ^ K {x 3 )X\x 2 ) 

-r K (l - r 2 2 )x 3 $ p K (x 3 )X L (x 2 ) + r K [2r\x 2 + (1 - r\)x 3 ) <$>° K {x 3 )X L {x 2 ) 
xE 6 (tf ) )hf ) (x 1 ,x 2 ,x 3 ,b 1 )} , (84) 

16nC F ^2N c J\dx] M^b^i A) 

x{[{l-2rl){l-x 2 )^ K ,{x 3 )X\x 2 )+C x l -r 2 2 ^ K ,{x 3 )X t {x 2 ) 

-r K .(l - r 2 2 )x 3 $ s Kt (x 3 )X L (x 2 ) + r K « (2r 2 2 (l - x 2 ) + (1 - 2r 2 2 )x 3 ) ¥ K ,(x 3 )X L (x 2 ) 

xE 4 {t { J ) )h i j\x 1 ,x 2 ,x 3 ,b 1 ) 

- \x 2 + (1 - 2r 2 2 )x 3 )<5> K *(x 3 )X L (x 2 ) + r 2 2 (2r K *¥ K ,(x 3 ) - U> K ,(x 3 ))X\x 2 ) 
-r K *(l - r 2 2 )x 3 ^(x 3 )X L (x 2 ) - r K * [2r\x 2 + (1 - 2r\)x 3 ) & K «{x 3 )X L {x 2 ) 
xE^tf^fix^x^xzM)} , (85) 



Cxl^C F ^2N c J^[dx} J™ b 1 db 1 ^ B (x 1 , 6i) 

x{[(l -x 2 + (1 - 2r 2 2 )x 3 )$^(a;3)X L (x 2 ) +Cxr 2 2 (2r^$^( a ;3) - ^ K *(x 3 ))X\x 2 ) 
-r K .{\ - r 2 2 )x 3 ^(x 3 )X L (x 2 ) - r K , (2r 2 (l - x 2 ) + (1 - 2r 2 )x 3 ) ¥ K ,(x 3 )X L (x 2 ) 
x^(4 1) )/i[ 1) (xi,i 2) x 3) 6i) 

- "(1 - 2r 2 )x 2 $ x *(x 3 )X L (x 2 ) + ir 2 2 $^(x 3 )X*(x 2 ) 

-r*.(l - r 2 ):r 3 $^(:r 3 )X L (:r 2 ) + rjr . (2r 2 x 2 + (1 - 2r 2 )x 3 ) & K «(x 3 )X L (x 2 ) 
xE 6 (t^)hf\x 1 ,x 2 ,x 3 ,b 1 )} , (86) 



16nC F] /2N c £[dx] J™ b 1 db 1 ^ B (x 1 , W] 
xr 2 { [(1 - r 2 )(l - :r 2 )^(:r 3 )X y (:r 2 ) 



Kx\r K *((l ■ riW K .[x 3 ) (1 r 2 )vI/;.. ! .r: i )).V'i.r 2 ) ^.i/;; )/// : (.r 



+ 



(1-^)X 2 ^.(X 3 )X^X 2 ) 
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M XK * 

JVI N6 



(1 - r 2 )*^) - \r K *{{\ + rlW K ,{x,) + (1 - r^fo)) X T (x 2 ) 
-2r K *((l + r 2 2 )x 2 + (1 - r 2 2 )x^ v K , (x 3 )X v (x 2 )} E A {tf)hf (x u x 2 , x 3 , h) } , (87) 

-( x l6nC F J2N c f\dx] H b 1 db 1 ^ B (x 1 ,b 1 ) 
v Jo Jo 

xr 2 { [(1 - r 2 2 )(l - x 2 )^ T K ,(x 3 )X v (x 2 ) 

~Cx ((1 - rl^Ux,) - \r K *{{\ + r 2 2 )$^(x 3 ) + (1 - r 2 )*^)) X r (* 2 ) 

-2r^((l + r 2 )(l -x 2 ) + (1 -rl)x 3 )^ v K .{x 3 )X v {x 2 )]E 6 {t^)h^{x 1 ,x 2 ,x 3 ,b 1 ) 
+ [(l-r 2 2 )x 2 V T K ,(x 3 )X v (x 2 ) 

+±r K *((l + r 2 2 )$ v K *(x 3 ) - (1 -r 2 )^(a:3))X T ( a : 2 )] J E 6 (i 2) )4 2) ( a ; 1 ,x 2 ,X3,6 1 )}88) 



32irC F ^2N c j\dx\ J™ M&i'M^iA) 

xr 2 { [(1 - ^^(^X^) - Cx^r x * (<t>^(rr 3 ) - (1 + 2r 2 2 )^(^)) X T (x 2 ) 
xE^t^h^ (x u x 2 , x^h) 



+ 



x 2 ^ K ,{x 3 )X v {x 2 ) - [& K .(x 3 ) - ^r K *($ v K *(x 3 ) + (1 + 2r 2 2 )^ K ,{x 3 )) )X T {x 2 ) 



1 



-2r x *((l + 2r 2 )x 2 + X3)$^(x 3 )X y (x 2 )] J B 4 (4 2) )4 2) (x 1 ,a; 2 ,X3,6 1 )} , 



(89) 



= -Cx^TrC F J2N c f\dx\rb 1 db 1 ^ B {x 1 M) 
v Jo Jo 

xr 2 {[(l-x 2 )^(x 3 )X v (x 2 ) 

-Cx (^(x s ) - ^T*.($k.(x 3 ) + ( X + 2r 2 2 )^(^3))) X T (x 2 ) 
-2r K .((l + 2r 2 )(l - x 2 ) + x 3 )$ a K .{x 3 )X v {x 2 )]E 6 (t$ ) )h$\x 1 ,x 2 ,x 3 ,b 1 ) 



+ 



~x 2 $ T K *(x 3 )X v (x 2 ) - -r K . (<t>^(x 3 ) - (1 + 2r 2 )^(:r 3 )) X T (x 2 ) 
xE 6 (4 2) )4 2) (a; 1 ,x 2 ,X3,6i)}, 
with = +1 an d Cxci = — 1 i n & U the above factorization formulas. 



(90) 



A.2 5 



(XcO^c)^ 



We present the spectator amplitudes for the B — > (x c0 , rf c )K^ decays below, where the symbol 
X represents the charmonium Xco or r\ c : 



Ml K = l6nC F ^m c j\dx] J™ hdb^Bix^h) 
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x{ [(1 - x 2 )$ K (x 3 )X v (x 2 ) - ( x r 2 2 (2r K <5> p K (x 3 ) - U K (x 3 ))X s (x 2 ) 

+r K (l - r 2 2 )x 3 ^ K (x 3 )X v (x 2 ) - r K (2r 2 2 (l - x 2 ) + (1 - r 2 2 )x 3 ) <S> p K (x 3 )X v (x 2 ) 

xE 4 (t { a ) )h { j\x 1 ,x 2: x 3: b 1 ) 

- \x 2 + (1 - 2r 2 2 )x 3 )<& K {x 3 )X v {x 2 ) + r 2 2 {2r K ^ K {x 3 ) - U> K (x 3 ))X s {x 2 ) 
-r K (l - r 2 2 )x 3 ^ K (x 3 )X v (x 2 ) - r K {2r\x 2 + (1 - r 2 2 )x 3 ) ^ K (x 3 )X v (x 2 ) 
xE^hfix^x^xM) , (91) 

(xl67cC F ^2N c £[dx] J™ hdh^BixiM) 

x{[(l - x 2 + (1 - 2rl)x 3 )^ K {x 3 )X v {x 2 ) - ( x r 2 2 (2r K <5> p K (x 3 ) - U K (x 3 ))X s (x 2 ) 
-r K (l - r 2 2 )x 3 ^ K (x 3 )X v (x 2 ) - r K (2r 2 2 (l - x 2 ) + (1 - r 2 2 )x 3 ) ^ K (x 3 )X v (x 2 ) 
x£ 6 (4 1) )/if ) (a; 1 ,a; 2 ,X3,6i) 

- ~x 2 *£> K (x 3 )X v (x 2 ) + r 2 2 (2r K <S> p K (x 3 ) - U K (x 3 ))X s (x 2 ) 

+r K (l - r 2 2 )x 3 ^ K (x 3 )X v (x 2 ) - r K [2r\x 2 + (1 - r\)x 3 ) <S> p K (x 3 )X v (x 2 ) 
xEs{tf)hf{ Xl ,x 2 ,x 3 M)}, (92) 

16nC F ^2N c J^[dx] jf°°l)A$B^i,&i) 

x{[(l -x 2 )$^(x 3 )X ?, ( a : 2 ) - ( x r 2 2 (2r K *<S> s K «(x 3 ) - ^ K ,(x 3 ))X s (x 2 ) 

+r K .(l - r 2 2 )x 3 ¥ K *{x 3 )X v {x 2 ) - r K * (2r 2 2 (l - x 2 ) + (1 - r 2 2 )x 3 ) $ s K *(x 3 )X v (x 2 ) 

xE 4 {t ( p)h i J\x 1 ,x 2 ,x 3 ,b 1 ) 

- \x 2 + (1 - 2r 2 2 )x 3 )$ K *(x 3 )X v (x 2 ) + r 2 2 {2r K *<$> s K {x 3 ) - h K *(x 3 ))X s (x 2 ) 
-r K .(l - r 2 2 )x 3 & K *{x 3 )X v {x 2 ) - r K , (2r 2 2 x 2 + (1 - r 2 2 )x 3 ) <Z> s K *(x 3 )X v (x 2 ) 
xE 4 {tf)hf{x u x 2 ,x 3 M)} , (93) 

C x 16nC F ^2N c j\dx\ J™ b 1 db 1 ^ B (x 1 , h) 

x{[(l -x 2 + (1 - 2r 2 2 ):r3)$^(:r3)X^ 2 ) - Cxr 2 (2r^$^ (x 3 ) - (x 3 ))X s ( a ; 2 ) 
-r*.(l - r 2 )^*^)^^) - r*. (2r 2 (l - x 2 ) + (1 - r 2 )x 3 ) & K .(x 3 )X v (x 2 ) 
xE 6 (t { J ) )h { J\x 1 ,x 2 ,x 3 ,b 1 ) 
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with £ 



XcO 



x 2 ^(x 3 )^(x 2 ) +r 2 2 (2r x *$^( a ;3) - -<$> K *{x 3 ))X s {x 2 ) 
+r K *(l - r 2 2 )x 3 <l> t K ,(x 3 )X v (x 2 ) - r K * (2r 2 2 x 2 + (1 - rj)x 3 ) $ s K *(x 3 )X v (x 2 ) 
xE 6 {t ( p)h < a ) {x l ,x 2 ,x 3 ,b l )} , 

+1 and ( Vc = — 1. 



(94) 
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